Competition between the Mott transition and the Anderson localization in ID 

disordered interacting electron systems 
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The competition between the Mott transition and the Anderson localization in one dimensional 
electron systems is studied based upon the bosonization and the renormalization group method. 
The beta function is calculated up to the second order in the strength of diagonal disorder by using 
a replica trick. It is found that the sufflciently strong forward scattering by random impurities 
destroys the Mott-Hubbard gap, and the backward scattering gives rise to the Anderson localization 
for the resulting gapless state. On the other hand, if the Umklapp interaction is strong enough, the 
Mott insulating state still overwhelms the Anderson localization. 



There has been continuous interest in the studies of 
correlated electron systems in the presence of disorder 
Although in the case of weak localization the ef- 
fect of electron-electron interaction has been extensively 
studied in the Hartree-Fock approximation, it is not well 
understood how the strongly correlated electron systems 
such as the Hubbard model are affected by the presence of 
disorder |^,^. In particular, the effect of disorder on elec- 
tron systems with the excitation gap, such as the Mott 
insulator, has received considerable attention recently. 

In one-dimensional(lD) electron systems, some elab- 
orate techniques like bosonization and renormalization 
group were successfully applied to disordered systems 
and revealed the role of electron-electron interaction 
for the localization-delocalization transition In 
the previous studies, however, the Umklapp interaction 
which is essential to cause the Mott transition was not 
taken into account. To investigate the disorder effect on 
the Mott insulating phase, it is crucial to study how the 
Umklapp interaction should affect low-energy properties. 

Motivated by this, we investigate ID disordered inter- 
acting electron systems putting emphasis on the compe- 
tition between the Mott transition and the Anderson lo- 
calization. We clarify the interplay between the disorder 
and the Umklapp interaction by using the bosonization 
and the renormalization group method. It is found that 
the Mott-Hubbard gap collapses and a gapless metallic 
state appears when the random forward scattering due 
to impurities is strong enough. The impurity backward 
scattering then drives this metallic state to the Ander- 
son localized state. We should mention here a recent 
numerical work for the ID Hubbard model with disor- 
der, which indicates the transition from the Mott insu- 
lator to the Anderson localization We note that our 
field-theoretic approach not only provides the knowledge 



complementary to the above numerical work, but also re- 
veals a new feature in the role played by the forward and 
backward impurity scatterings. 

We start by introducing the effective Hamiltonian to 
describe ID interacting electron systems with random po- 
tentials. Applying standard bosonization rules, we write 
down the effective Hamiltonian in field theory limit [0 , 
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Hdis = \ - dxr]{x)dx(t>c{x) 



dx{C(x)e^(^'^<=(^)+2'=''^Hr(.9(a;)) + h.c.}. (4) 



Here 0c is a boson phase field for the charge degrees of 
freedom and He is its canonical conjugate momentum 
field, and S = Akp — 27r. The [/-term in eq.(^) rep- 
resents the Umklapp interaction, which may cause the 
Mott transition. For the spin degrees of freedom, we 
have used non-abelian bosonization |^,^ to preserve 
SU(2) symmetry: J l(R) is the left(right)-going current 
operator of level-I SU(2) Kac-Moody algebra, and g{x) 
is a fundamental representation of SU(2) Lie algebra. 
The A-term in eq.(|^) is a marginally irrelevant inter- 
action which arises from SU(2)xSU(2) symmetry. We 
have introduced real and complex random fields ?7(x) and 
^{x) for forward and backward scatterings by impurities. 
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respectively, which obey the gaussian distribution law, 
ivixUx')) = D,S{x - x'), {Cix)ax')) = D^5{x - x'), 
and {ri{x)) = {£,{x)) = 0. In addition to the above inter- 
actions, we should consider the following term 



dxA(x)dx(- 



(5) 



where dxdc = He and A{x) is a random gauge field with 
the gaussian distribution, {A{x)A{x')) = Da6{x — x'), 
{A{x)) — 0. As will be seen later, this term is generated 
in the process of renormalization due to the backward 
scattering by random potentials. We note that apart 
from the [/-term, the charge part of the Hamiltonian 
is similar to that of the 2D random phase sine-Gordon 
model for which randomness is introduced for the 

2D plane in contrast to the present ID system. 

We now consider the quenched disorder with the use 
of a replica trick. Introducing n species of replicas and 
integrating over the random variables [p^ , we have the 
following effective action arising from disorder, 

Sms^-'^ j j dT J dT'Y,9x(t'lix,T)dx^,{x,T') 

j dx j dT j dT'Y,^^drCbl{x,T)dr<l)i{x,T') 
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[dxfdrf dT'J2ir{g\x,T))tr{g'{x,r')) 



X COS V27r((/>*(a;,r) - 0^(x,t')), 



(6) 



where i, j are replica indices. In the following argument, 
we calculate the beta functions up to the second order 
in D^^ri,A- To this end, it is quite useful to exploit the 
operator product expansions for the U(l) gaussian model 
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with z = X + iv^T, z — X — iv^T, and also those for the 
level-1 SU(2) Wess-Zumino-Witten model fT^-jl^, 

Jlix, r) Ji(0, 0) ~ % + — J£(0, 0) + • • ., (9) 
2:^ z 

J£(x,T)g(0,0)~^.g(0,0) + ---, (10) 
tr(5(a;, r))tr(5(0, 0)) ^ |z| Jl(0, 0) • Jr{0, 0) + ^ 



^(z2Jl(0,0)- Jl(0,0) 
z^Jh{Q,0)-Jh{0,0)) + - 



(11) 



with z = X + ivsT, z = X — ivsT, and t"", the generator 
of the SU(2) Lie algebra. By expanding the action in 
terms of the interactions and using the above operator 
product expansions, we end up with the following scal- 
ing equations up to the second order in the strength of 
randomness and the lowest order in U and A after taking 
the replica limit n —> 0, 



dD. 



dl 

dn 



V 



dl 

dPA 

dl 

dU 
H 
d~X 
'dl 
dKc 
dl 
dvc 
~dl 
du 

U 



(2 - ife - 3A)D5, 



A2 ~ 

De, 

2 ^' 

-27rii:2i7Vo(fc) 
-hi ' 



2u 



(12) 
(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 



where b^^r,,A = D^j^^a/vI, U = U/vc, A = A/w^, 



9{u) 
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-00 y/l + u^y^l+y^)^^/^ 



(20) 



and Jo{x) is the Bessel function. We see from eq.(|l4[) that 
DA-term is indeed generated by the second order contri- 
bution of ZJj-term, even when Da is initially equal to zero 
as mentioned before. We can see that Da does not cou- 
ple to the Umklapp term, so that the random gauge field 
may not cause any essential change in the Mott insulat- 
ing phase. The effect of the random gauge field manifests 
itself in the metallic phase away from half-filling, as will 
be discussed later again. Up to the first order in D^ and 
the zeroth order in Z?,, and Da, these scaling equations 
coincide with those obtained by Giamarchi and Schulz 
except that they did not take into account the Umklapp 
interaction JTo| , pO[ . 

In what follows, we restrict our arguments to the case 
of half-filling. We then put (5 = and Jo(0) = 1 in 
eq.(^). In the absence of randomness, the Umklapp in- 
teraction is relevant, resulting in the low-energy fixed 
point of the Mott insulator with the charge excitation 
gap. Let us first study the effects of the random for- 
ward scattering D,j on this Mott insulating phase. For 



2 



a while we neglect the backward scattering, by putting 
= and Da = (I)^-term is generated by ij^-term) . 



Then A in is scaled to zero after the renormaliza- 
tion, preserving spin excitations still massless (the inter- 
play between A and D{ will be discussed later). Also, 
from eq.(|l^), we have Dr]{l) — i),,(0)e'. Substituting 



this expression into eq.(|l5[), and solving eqs.(|l5[) and ( |17|) 
numerically, we obtain the renormalization flow for U 
and Dr, as shown in Fig. 1. We find that for suffi- 
ciently large values of £',,(0), the Umklapp term, IJ{1), 
becomes irrelevant and thus the Mott-Hubbard gap dis- 
appears, resulting in massless charge excitations. Note, 
however, that this massless state is not a conventional 
Tomonaga-Luttinger liquid but is a disordered metallic 
state for which some spatial correlation functions show 
exponential decay, as discussed in On the other 

hand, if U{Q) is sufficiently large compared to Djj{0), we 
can see that U{1) is scaled to a strong-coupling value as 
I — > +00, and thus the Mott-Hubbard gap should still 
persist. The transition between the Mott insulator and 
the gapless metallic phase is of Kosterlitz-Thouless type 
Also, the numerical results indicate that at the 
critical point the magnitude of the Mott-Hubbard gap 
is roughly proportional to the inverse of the correlation 
length, l/ln(l)-i(0)). 

Having noticed that the strong forward scattering can 
drive the Mott insulator to a disordered metallic state, 
let us now discuss the effects of the backward scattering 
which may bring about the Anderson localization. 
The scaling equation (|l^) implies that for the repulsive 
U {Kc < 1), is renormalized to a larger value even if 
its initial value is small, and hence the backward scatter- 
ing becomes relevant. Therefore, if there is no intermedi- 
ate fixed point between the weak-coupling and strong- 
coupling regimes, the low-energy fixed point is classi- 
fied by C7* = and i)| +oo, or by [7* ^ oo and 
_D| +00. In both cases, the fixed-point is identified 
with the insulator, since the Drude weight D = v^Kc is 
scaled to zero by U oi D^. In order to see whether this 
is the Mott insulating state or the Anderson localized 
state, we examine the behavior of the compressibility. 
The compressibility is given by Kcjvc-, which satisfies, 



dl \ Vr 



(21) 



as easily seen from eqs. 

(0) and (IH). Thus if U is irrel- 
evant, the compressibility takes a finite value, resulting 
in the fixed point of the Anderson localization. On the 
other hand, if U is relevant, the compressibility vanishes 
at the low-energy fixed point, characterizing the Mott in- 
sulator. Therefore when the random forward scattering 
Dri{Q) is in the region where [/ — > in Fig. 1, the back- 
ward scattering -0^(0) drives the system to the Ander- 
son localized state, whereas for sufficiently large values 
of U{0) compared to D,j{0), the Mott insulating state still 



persists. Therefore, we end up with the conclusion that 
there can occur the transition driven by disorder from 
the Mott insulator to the Anderson insulator for ID dis- 
ordered electron systems, although the Mott insulator 
still overwhelms the Anderson localization if the Umk- 
lapp interaction U(0) is large enough. This is consistent 
with the numerical results for the disordered Hubbard 
model |]. 

The above characteristic behavior found for correlated 
electron systems shows sharp contrast to the results for 
ID spinless fermion systems A remarkable point 

for the spinless fermion model is that the Mott transi- 
tion is always accompanied by the order of 2fci^^-charge 
density wave (CDW). Since this state can be mapped to 
the Ising ordered state, one can apply Imry-Ma's state- 
ment |2^ that infinitesimally small disorder destroys this 
Ising ordered state. This indeed leads to the conclusion 
that the Mott insulator for the spinless fermion model is 
changed to the Anderson localized state even if infinites- 
imally small disorder is introduced Obviously, this 
argument cannot be directly applied to the Mott insula- 
tor for the present electron systems, which is not accom- 
panied by the CDW order. 

The difference between the spinless model and the elec- 
tron model can also be seen in the scaling equations 
(|l3| ) and (|l^): for the spinless model g{u) is replaced by 
{T{Kc - l/2)/T{Kc))'^, where T{x) is the gamma func- 
tion. While g{u) is always finite for the electron model, it 
diverges for Kc < 1/2 for the spinless fermion model, and 
the renormalization group procedure breaks down. Thus 
as Kc approaches this singular point, £)„ becomes large, 
and at last the right-hand side of cq.(|l5|) becomes neg- 
ative, even for infinitesimally small D^^jj{Q). Note that 
the Mott transition for the pure spinless model occurs at 
Kc = 1/2. As a result, the Umklapp term U becomes 
always irrelevant and the Mott gap closes for a disor- 
dered spinless model. The above comparison with the 
spinless fermion model naturally leads us to claim that 
the stability of the Mott insulator for electron systems 
against weak diagonal disorder may be closely related to 
the absence of the long-range order. 

So far, we have been concerned with the charge degrees 
of freedom. Here we briefly mention the effect of disor- 
der on the spin degrees of freedom. The scaling equa- 
tion ([l^ ) indicates that the backward scattering may 
renormalize A to a negative large value, hence freezing 
the spin degrees of freedom. This is essentially the same 
as that observed for the case away from half filling . 
Therefore, the Mott-insulating phase with frozen spins 
may be realized if the spin degrees of freedom is frozen 
after the charge-gap formation due to the Umklapp term 
(the large-?/ region in Fig.l). Also, for the small-C/ re- 
gion, the Anderson localized state with frozen spins may 
be expected. 

Finally some comments are in order for the random 
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gauge field A{x) generated by the impurity backward 
scattering. Although this term does not cause a seri- 
ous change in the half-filled Mott insulator, it plays an 
important role in a delocalized phase. To see this clearly, 
we consider the system away fi'om half-filling. We can 
then put [/ = in the scaling equations because the 
Umklapp interaction is irrelevant. We see from eq.(p^) 
that if the electron-electron interaction is strongly attrac- 
tive and the condition, 2 — Kc — 3A < 0, is satisfied, the 
backward scattering becomes irrelevant and then a 
delocalized metallic phase realizes. It has been naively 
expected that in this phase the correlation for singlet 
superconductivity is dominant [@-|o|. It turns out, how- 
ever, the random gauge field may change properties of 
this phase. The random gauge field A{x) as well as ijix) 
can be incorporated into the shift of the phase fields . 
Defining new phase fields as, 

ec{x) = ec{x)+A{x), (22) 
A{x)^-^J^ r dx'A{x'), (23) 

etc., we can cast the charge part of the Hamiltonian into a 
gaussian type. Although the random field, A{x)^ is now 
eliminated from the Hamiltonian, it may change long- 
range behaviors of some correlation functions. For ex- 
ample, the random gauge field A{x) brings about the 
exponential decay of the correlation function for singlet 
superconducting pairing: 

(ctl(x)cx«(x)c|^(0)c|^(0)) ^ —L^\^-^'^\x\e-'Mf. 

\X\ -^<: 

(24) 

This is due to the random phase shift caused by the scat- 
tering due to random gauge fields. Eq.(|2^) implies that 
the susceptibility for singlet superconductivity does not 
diverge in the thermodynamic limit, Xss{k = 0,0;) 
uj^/^" . Thus the fiuctuation toward superconductivity is 
much suppressed by random gauge fields. 

In conclusion, the transition between the Mott insu- 
lator and Anderson insulator occurs according to the 
strength of disorder and the Umklapp interaction. We 
have found that the random forward scattering by impu- 
rities plays a central role to change the Mott insulating 
state to a metallic state, whereas the backward scatter- 
ing drives the resulting metallic state to the Anderson 
localized state. 
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FIG 1. The renormalization flow for U and Z?.,,. The 
flow for some initial values of U is shown: Kc{Q) = 0.8, 
D,,(0) = 0.05. The transition point exists around U{Q) — 
0.0463. 
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